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Compared with direct-gap semiconductors, the valley degeneracy of silicon and germanium opens 
up new channels for spin relaxation that counteract the spin degeneracy of the inversion-symmetric 
system. Here the symmetries of the electron-phonon interaction for silicon and germanium are iden- 
tified and the resulting spin lifetimes are calculated. Room-temperature spin lifetimes of electrons 
in silicon are found to be comparable to those in gallium arsenide, however, the spin lifetimes in 
silicon or germanium can be tuned by reducing the valley degeneracy through strain or quantum 
confinement. The tunable range is limited to slightly over an order of magnitude by intravalley 
processes. 

PACS numbers: 71.70.Fk, 72.25.Dc, 72.25.Rb, 72.10.Di 



I. INTRODUCTION 



The favorable material properties of silicon have per- 
mitted it to dominate the microelectronics industry for 
over half a century, however a new genre of spintronic 
semiconductor devices^— in which spins of electronic 
carriers are manipulated instead of a charge current, re- 
quires long spin transport lengths and coherence times. 
Although spin injection into nonmagnetic semiconduc- 
tors was demonstrated over a decade ago,— ~— the recent 
success at injecting spin-polarized current into silicon^— 
suggests incorporation of semiconductor spintronic de- 
vice concepts into hybrid silicon device architectures. Po- 
larized spins relax in semiconductors because the spin- 
orbit interaction entangles orbital and spin degrees of 
freedom, and thus ordinary scattering from defects or 
lattice vibrations leads to a loss of spin coherence. In 
materials without inversion asymmetry the entanglement 
of spin and orbit manifests as an effective momentum- 
dependent (internal) magnetic held, causing spin pre- 
cession and D'yakonov-Perel' spin relaxation. Within 
inversion-symmetric materials, such as silicon, the in- 
ternal magnetic field vanishes, but scattering between 
states with spin-orbit entangled wave functions leads to 
Elliott- Yafet spin relaxation. The spin coherence times 
in silicon are long at low temperature, and the spin- 
orbit interaction and lattice symmetry reduces spin re- 
laxation rates relative to optically accessible (direct-gap) 
semiconductors The silicon band structure, how- 
ever, has multiple valleys that permits low-energy scat- 
tering of electrons by large momenta, which allows the 
Elliott- Yafet process to be more effectived Numerical 
calculations that include these effects have been success- 
ful at explaining the spin lifetime in silicon as a function 
of temperature ! ' 18 Tuning the spin lifetime in inversion- 
asymmetric semiconductors with a single direct gap have 
largely focused on the electric-field-induced Rashba spin- 
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orbit interaction that shortens the spin lifetime 
addressed is the potential for new methods of tuning the 
spin lifetime associated with the valley degeneracy of the 
semiconductor . 

Here we trace that origin of the intrinsic spin lifetime in 
silicon to the spin flips associated with large momentum 
transfer events, and disentangle the intervalley contri- 
bution to spin relaxation from the intravalley contribu- 
tion. As the different processes involve different momen- 
tum transfers, and for the intrinsic spin-relaxation rate 
the source of that momentum transfer will be electron- 
phonon scattering, the various contributions will be as- 
sociated with specific regions of the phonon dispersion 
curves. Due to the high symmetry of the crystal lat- 
tice, several processes that might have been expected to 
contribute will be forbidden by symmetry. We provide 
a full symmetry analysis of the various contributions to 
the spin-relaxation rate from phonon-mediated scatter- 
ing. The separation of spin-relaxation mechanisms by 
momentum transfer also permits a direct calculation of 
the tuning of spin lifetime possible by splitting the en- 
ergies of the electron valleys and thus suppressing some 
of the intervalley scattering effects. Reducing the val- 
ley degeneracy in silicon, through applied strain or the 
growth of pseudomorphic SiGe quantum wells, will re- 
duce the effect of the spin-orbit interaction on electron 
scattering, lengthening the spin coherence time and spin 
transport length. We find an effective tuning range of 
approximately one order of magnitude. 



II. ELECTRON-PHONON SCATTERING 

The intrinsic spin-relaxation time is determined by 
electron-phonon scattering. For one-phonon absorption 
(+) and emission (— ) processes, the scattering probabil- 
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FIG. 1. (Color online) Electron spin-relaxation time (a) and 
mobility (b) in bulk silicon shown by black solid lines. In- 
dividual scattering processes are intravalley acoustic process 
(dashed line), intravalley optical (rj 5 ) process (plus symbol), 
intervalley A processes (closed symbols), and intervalley E 
processes (open symbols). 



FIG. 2. (Color online) Electron spin-relaxation time (a) and 
mobility (b) in bulk germanium shown by black solid lines. 
The dashed line shows the intravalley acoustic phonon con- 
tribution and the plus symbol shows the intravalley optical 
phonon contribution. The open symbols show the intervalley 
X processes for three effective phonon energies. 



ity from k to k' is 

M± CT (k',k) = A/W,n qT l 



(1) 



where H 1 ^ is the time-independent part of the electron- 
phonon interaction Hamiltonian corresponding to ab- 
sorption or emission, q = k' — k, a labels the spin 
state, and n q is the phonon occupation number. To 
evaluate CT (k',k) for various types of scattering pro- 
cesses in the material, we use a first-nearest-neighbor 
sp 3 tight-binding model (TBM) with on-site spin-orbit 
interactions 2 ^* to obtain the wave functions, 
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where N is the number of unit cells, j labels unit cells, a 
labels the two basis atoms within a unit cell, I labels the 
atomic orbital bases, x is a two-component spinor, s is 
the spin index, and R JQ is the position vector of atoms. 
We choose the spin quantization axis to be aligned with 
the z axis and determine the coefficients c a ; s by maxi- 
mizing the expectation value of the spin operator (S z ). 

In the spherical band approximation, we express the 
matrix elements in terms of the deformation potentials 
up to the first order in 5q = (k' — k/) — (k — k^), 22 
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where p is the density, V is the crystal volume, lo is 
the phonon frequency, and k^ and k^ are the wave 
numbers of the initial and final valley minima, respec- 
tively. For Si, there are six valleys on the A axes, e.g., 
kg = (27r/a)(0.85, 0, 0), and for Ge, there are four valleys 



at the L points, e.g., (7r/a)(l, 1, 1). A spherical averag- 
ing around the valley minimum is carried out for eval- 
uating D\. We have assumed only the following types 
of matrix elements are nonzero, (</>as(R.)|f§|<^ap 

(^(R')llf I'M*)), and (^'(ROIa&l&fCH-)). and 
all have the same magnitude. Since the atomic potential 
is not explicitly known in the TBM, we can only deter- 
mine the relative strengths of the processes. The overall 
magnitude is later determined by fixing the mobility to be 
1450 cm 2 /Vs in Si and 3800 cm 2 /Vs in GeS 3 . Calculated 
deformation potentials for different phonon processes are 
summarized in Table HI The relative strengths of different 
processes are consistent with the semiempirical values in 
the literatures 2 .! 2 ^ 

Various electron-phonon scattering rates are computed 
as follows. For the intravalley acoustic process, a linear 
phonon dispersion, uj — c|q| is used and the scattering 
rate including both absorption and emission is 
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where m = (totto 2 ^) 1 / 3 is the averaged effective mass, 
and c = 3/(l/c| + 2/c 2 ^) is the averaged speed of 
sound. We use p = 2329 kg/m 3 , iul = 0.9163m e , 
tot = 0.1905to £ , cl = 8500 m/s, and ct = 5900 m/s for 
Si, and p = 5323 kg/m 3 , = 1.59TO e , to^ = 0.0823m e , 
c L = 4900 m/s, and c T = 3500 m/s for GeS 3 - A ther- 
mal averaging of the initial electron energy is carried out 
with the Boltzmann distribution, appropriate for non- 
degenerate systems, 

(g(E)) T ee ^=±_ j( g{E )E^e- E ' T dE. (5) 

In this regime the spin lifetime is independent of the car- 
rier density; for degenerate carrier densities the spin life- 
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TABLE I. Deformation potentials. The unit for Do's is 
eV/A, and for LVs is eV. For intravalley acoustic processes, 
D A = 3.1 eV, D A ,n = 0.0016 eV for Si, and D A = 3.8 eV, 
D A ,ti = 0.032 eV for Ge. For spin-flip processes, the super- 
script xy indicates that both the initial and final valleys are in 
the x-y plane, and the superscript z indicates that the initial 
and final valleys are separated along the z direction. T UQ is 
the effective phonon frequency in Kelvin. The deformation 
potentials listed here do not include the valley degeneracy of 
final states. The parentheses indicate that the D\ term is 
negative. 

Phonon T„ (K) D D, D£ n D^ n Z3g, n Df iU 
T+ 730 0.34 0.14 0.19 
A2+A5 700 4.5 (3.2) 0.031 0.01 
Ai 210 0.04 0.028 0.04 
1 A 5 140 2.7 0.009 0.0015 
E1 + E2 630 3.2 (2.3) 0.03 0.11 0.044 0.17 
Ei + E 3 500 3 (1.1) 0.018 0.075 0.026 0.097 
S 3 + E 4 210 0.0083 2.2 0.0083 0.041 0.0059 0.058 
430 3.5 4.6 
Ge X 4 390 0.24 1 0.24 0.16 0.17 
Xi 340 3.8 3.4 0.08 0.054 0.12 0.089 
X 3 120 2.7 0.11 0.1 



TABLE II. Selection rules with the inclusion of time-reversal 
symmetry without spin. For electron representations, Ait is 
Ai at ko transformed to a valley on a perpendicular axis to 
ko, and Lit is L\ transformed from one L valley to a dif- 
ferent valley. The phonon representation Ai(2ko) becomes 
A 2 (2ko — (47r/a)(l, 0, 0)) in the reduced Brillouin zone scheme 
used in Table [I] 
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TABLE III. Selection rules with the inclusion of time-reversal 
symmetry with spin. The spin-flip part of Ai (Aj in Table [J) 
in Si and of F^ 5 in Ge are also forbidden by time reversal. 
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time will be shorter than for nondegenerate carrier den- 
sities as the distribution function spreads out to a larger 
k range for which the spin mixing upon scattering will 
be larger. For the intravalley optical process and each 
intervalley process, we use an effective phonon frequency 
Wo , listed in Table Q] as temperature T wo i 24 i 25 The zeroth- 
order and the first-order scattering rates are 
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p ( ±> = y/E ± hu 6{E ± TUjJq), 

p { ±> = (2E ± Tiuj )\/E ± Tiujq 6(E± hw 



I) 



(8) 
(9) 



The final electron spin-relaxation times including all scat- 
tering processes for Si and Ge are plotted, respectively, 
in Figs. [1] and [5] as a function of temperature. 



III. DISCUSSION 



representation at the conduction-band minima in Si (Ge) 
becomes Ag (Lj) instead of Ai (L^")^l We analyzed the 
selection rules including time-reversal symmetry using 
the same subgroup technique developed in Refs. [2fj| and 
[H and list the results in Table Mil The allowed phonon 
representations at q (right-hand side of the equations) 
are obtained from the characters at k and k' (the two 
representations on the left-hand side of the equations). 
Details of our calculations of the selection rules are pre- 
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Our results can be qualitatively understood from the 
selection rules, derived from symmetry, that apply to 
scattering processes between the valley minima, i.e., 
Do's. The selection rules without spin-orbit coupling 
were discussed by Lax and Hopfield— and are listed in 
Table [II] When electron spin is included, the irreducible 



FIG. 3. (Color online) Spin-relaxation time at T = 300 K as 
a function of valley energy shift, AEo, shown by black solid 
lines. The energy shift is the energy offset of four (three) 
valleys relative to the remaining two (one) in Si (Ge). The 
symbols in panel (a) have the same meaning as in Fig. [1] and 
in panel (b) are the same as in Fig. [21 



4 



sented in the Appendix. 

In Si the spin-orbit interaction mixes spin up and down 
by about 1% in wave function and gives a spin- flip prob- 
ability about 10~ 4 . It can be seen from Fig. [T]that the 
£ phonons (also known as the / processes) dominate 
the mobility and the spin flip near room temperature. 
For mobility, A phonons (also known as the g processes) 
also contribute substantially, consistent with the selec- 
tion rules of A' 2 and £1 in Table [TT1 For spin flip, the £2 
and £3 phonons become allowed with nonzero spin-orbit 
interaction in addition to £1 (Table Ull]) . The intravalley 
optical (rj;) phonons and the A phonons remain for- 
bidden for spin flip by time reversal and, therefore, are 
not as effective as the £ phonons. Although the cou- 
pling strength for the low-energy £ phonon is weaker, 
this is somewhat compensated by the temperature de- 
pendence of the phonon distribution and it ends up that 
all £ phonons contribute approximately the same to spin 
flip near room temperature. 

In Ge the spin-flip probability due to spin-orbit inter- 
action is about one order of magnitude larger, but the 
number of possible phonon processes is reduced. So the 
spin-relaxation time is only about one order of magni- 
tude smaller than in Si as shown in Fig. [2] For mobil- 
ity, the intervalley X± and the intravalley optical (rj 5 ) 
phonons are more important due to the selection rules 
in Table [TTJ For spin flip, the high-energy (A4) phonon 
is allowed with finite spin-orbit interaction, but the low- 
energy (A 3 ) phonon is still forbidden by time reversal. 
Although the coupling to T^ 5 phonons is as strong as the 
coupling to Ai phonons, the spin-flip part is forbidden 
by time reversal and ineffective (Table ITTT|) . 

Now that the structure and symmetry of the spin- 
relaxation mechanisms has been clarified, the analysis of 
the effect of strain is straightforward. Strain (or quan- 
tum confinement) breaks valley degeneracy and can elim- 
inate multivalley scattering processes. In Si, a [001] 
strain can change the lowest-energy valley degeneracy 
from six to two (located on the same axis). In Ge, a 
[111] strain can yield just a single nondegenerate valley 
at the conduction-band edge. We have performed a sim- 
ple estimate that only takes into account the valley en- 
ergy shift, AEo, by modifying ±Tkjj —> ±hcu — AE and 
the initial electron distribution in each valley. As shown 
in Fig. |3l the spin lifetime averaged over all valleys can 
be lengthened substantially; 1% strain gives about 0.09 
eV shift in Si, and 0.16 eV shift in Gei^ Positive energy 
shift (AEo > 0) corresponds to the configuration that 
four valleys out of six arc shifted to higher energy in the 
case of Si, and three out of four in the case of Ge. Neg- 
ative shift reverses the ordering. For positive shifts, the 
intervalley £ processes in Si or the A processes in Ge 
can be completely eliminated. This tuning is eventually 
limited by the intervalley A and the intravalley optical 
(T25) processes in Si, and by the intravalley acoustic pro- 
cesses in Ge. For negative shifts, the elimination of the 
intervalley processes is only partial, so the tuning range 
is much smaller. 



IV. CONCLUSIONS 

We have presented a thorough symmetry analysis of 
the electron spin-phonon interaction processes for silicon 
and germanium, finding a spin lifetime for nondegener- 
ate carriers at room temperature comparable to those 
in III-V semiconductors when the scattering determin- 
ing the carrier mobility is dominated by phonons. How- 
ever, strain or quantum confinement can lift the valley 
degeneracy, which lengthens the spin lifetime substan- 
tially (over an order of magnitude at room temperature) . 
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Appendix A: Selection Rules 

In this Appendix we show the calculations of the char- 
acters of the product representations for obtaining the 
selection rules in Tables Ull and IIIII Each electron-phonon 
scattering process involves three subgroups of the wave 
vectors, k, k' and q. Instead of using the intersection 
group formed by the elements common to all three sub- 
groups and creating the corresponding character tables, 
the selection rules are derived using the existing charac- 
ter tables of the subgroups at these wave vectors. The 
two multiplication rules, with and without time-reversal 
symmetry, will be presented! 26 ' 28 We choose to compute 
the characters of the phonon representations at q from 
the products of the electron representations at k and k'. 
Once the characters at q are obtained, the decomposition 
into irreducible representations is done in the usual way 
using the subgroup of the wave vector q. 

TABLE IV. Double group characters at A = k (l, 0, 0). 



Class 


Ai 


A 2 


A 2 


Ai 


A 5 


A 6 


A 7 


(E\0) 


1 


1 


1 


1 


2 


2 


2 


(E\0) 


1 


1 


1 


1 


2 


-2 


-2 


(cf|o),(cf|o) 


1 


1 


1 


1 


-2 








2(C 4 |r) 


A 


— A 


-A 


A 





v^A 


-v^A 


2(C 4 \r) 


A 


—A 


-A 


A 





-v^A 


v^A 


2(iCt\r),2(iCl\r) 


A 


A 


-A 


-A 











2(iC 2 |0),2(iC , 2 |0) 


1 


-1 


1 


-1 












A = e ikoa/4 , E is 2tt rotation, iV A star a(C) = 1, 
A 6 = Ai ® D 1/2 , A 7 = A 2 ® D 1/2 

We first carry out calculations without time-reversal 
symmetry. To obtain the product character at q for each 
class C using the characters at k and at k', the usual 
character multiplication rules arc modified as follows. 



xf J (C) 



xKQ xl(C) 



N, 



q star k 



(C), (Al) 
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TABLE V. Double group characters at L = (■k/o)(1, 1, 1). 



Class 


, _|_ 

if 


if 


if 


if 


if 


_ _|_ 

if 


(£|0) 


1 


1 


2 


1 


1 


2 


(E\0) 


1 


1 


2 


— 1 


— 1 


— 2 


3(C 2 |r) 


1 


-1 







— i 





3(C 2 |r) 


1 


-1 





— i 


i 





2(C a |0) 


1 


1 


-1 


-1 


-1 


1 


2(C 3 jO) 


1 


1 


-1 


1 


1 


-1 


(i|T)Z 




±x(z) 






±x{z) 





Z can be any of the six classes shown above, 
if = Lf ® D1/2 =Lf® D 1/2 



TABLE VIII. Group characters at A = (27r/a)(l, 0, 0). 



Class 


Xj 


A2 A3 A4 


Ax star L 


(E\0) 


2 


2 2 2 


4 


(C||0) 


2 


2 -2 -2 





{C2\T),{C2\T + t X y) 





—2 2 


2 


2{iC 2 \Q) 


2 


-2 


2 


(E\t X y) 


-2 


-2 -2 -2 


4 


(C 4 \t x y) 


-2 


-2 2 2 





(C 2 \r + t xy ),(C' 2 \T) 





2-2 


2 


2(iC 2 \t xy ) 


-2 


2 


2 



= (a/2)(l,l,0) 



where i and j are the irreducible representations, and 
Nq s tar k(C) is the number of wave vectors that are un- 
changed by the class C, out of a set of nonequivalent k 
points. This set of nonequivalent k points is generated 
by the subgroup of q and is called "q star of k." The 
product character is simply zero if C is not a common 
class of all three subgroups. 



TABLE VI. Group characters at T = (0,0,0). 



Class 


r ± r ± 
1 1 1 2 


1 12 


1 15 


1 25 


Ar star A 


Ar star L 


(E\0) 


1 1 


2 


3 


3 


6 


4 


3(C 4 2 |0) 


1 1 


2 


-1 


-1 


2 





6(C 4 |t) 


1 -1 





1 


-1 


2 





6(C 2 |r) 


1 -1 





-1 


1 





2 


8(C 3 |0) 


1 1 


-1 











1 


(i\r) 


±1 ±1 


±2 


±3 


±3 





4 


Z(iCl\r) 


±1 ±1 


±2 


Tl 


Tl 


4 





6(iC 4 |0) 


±1 T l 





±1 


Tl 








6(iC 2 \0) 


±1 T l 





Tl 


±1 


2 


2 


8(«7 3 |t) 


±1 ±1 













1 



TABLE VII. Group characters at E = k (l, 1, 0). 



Class 


Ei 


E 2 


E3 E4 


As star A 


(E\0) 


1 


1 


1 1 


2 


(C 2 |r) 


A 2 


A 2 


-A 2 -A 2 





(iCl\r) 


A 2 


—A 2 


-A 2 A 2 


2 


(iC 2 |0) 


1 


-1 


1 -1 









A 2 = 


ik a/2 





Listed here are the character tables for the double 
groups representing the symmetry of the electron states 
at A for Si (Table EJ and at L for Ge (Table E}, and the 
tables for the groups at A, T, E, and X for phonons (Ta- 
bles |IV] and EEEini) • The numbers of invariant q-star-k 
points, Aq s tar k(C), are listed in the tables for phonons. 
Note that the group at X has 14 irreducible representa- 
tions, 14 classes, and 32 elements. Only the 4 physically 
admissible representations and 8 classes are listed in the 



table. The other 6 classes are not relevant because they 
have zero character for the 4 physical representations. 
Further simplification is possible because one can work 
with just 4 classes; the other 4 classes, corresponding 
to an additional lattice translation, t xy = (o/2)(l, 1, 0), 
have exactly the opposite characters. 

The relevant characters of the product representations 
are listed in Tables EHXI] for Si and in Tables EnTXini 
for Ge. Because the space group for the diamond struc- 
ture is non-symmorphic, some of the group elements con- 
tain a sublattice translation, r = (a/4) (1,1,1). Note 
that the complex phase factors (A) due to this transla- 
tion in some characters should always cancel out in the 
final character product (k' = k + q) and not affect the 
selection rules. We will ignore these explicit phase fac- 
tors in our character product tables. However, care needs 
to be taken with the product classes of L and L t when 
comparing to the classes at X. The product of two group 
elements with r can result in a lattice translation, t xy , 
which is explicitly included in the subgroup of X, but not 
in the subgroup of L. All physical representations at X 
are odd under t xy . The characters at L = (it /a)(\, 1, 1) 
are odd under t xy , but are even at L t = {■k/o){— 1, 1, 1). 
Therefore, the product class for (C2|r) at L and at L t is 
actually (C2\t T t xy ) at X. This is why the product of 
L\ and L\ t contains X3, but not 

Time reversal can add additional constraints to the 
selection rules if there exists a group element that con- 
nects the time-reversed process (— k' — > — k) to the orig- 
inal process (k — > k'). That is, we need an element Q 
from the subgroup of q that interchanges k and — k'. To 
incorporate the time reversal symmetry, the character 
multiplication rule is modified to the symmetric or anti- 
symmetric combination of the characters of the original 
process and the process connected via Q, 



XT±(C) = ^ XT( C ) ± XiiC^Nq star k (QC) (A2) 



where C 2 is the class of the square of elements in C and 
Q is the element that interchanges k and — k'. The pos- 
itive (symmetric) sign is used in the cases without spin 
and the negative (antisymmetric) sign is used with spin. 
Again, all the relevant characters for the product repre- 
sentations are listed in Tables HXI IXIIll The element Q 
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can be identified as the first entry in the column of QC 
when C is the identity class. For A phonons in Si and Y 
phonons in Ge, Q is simply the identity element. Upon 
closer examination, we also found that the spin-flip part 
of the Ai phonon process in Si and of the Y^ 5 phonon 
process in Ge are forbidden because the final state is ex- 
actly the time reverse of the initial state (e.g., \a) = |k t) 
and \ fa) = | — k J,)), 

(fa\H ep \a) = (falfHlf-^a) , (A3) 



where T is the time-reversal operator, T 2 = — 1 in the 
presence of half-integer spin, and H cp = TH^T^ 1 is 
the time-reversal invariant electron-phonon interaction 
Hamiltonian. 
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TABLE IX. Characters of the product representations at T. 



G 


QC 


N(QC) 


G 


x kn x (c ) 


Ai®i 

Xr 


\l Ai®Ai 

Xrj 


G 


X kn 6 (C ) Xr 


6 ®A 


6 A 6 (gA 6 

Xr- 


(E\0) 


(i\r) 





(E\0) 


1 


6 


3 


(E\0) 


2 


24 


12 


(C||0) 


(iCj\r) 


4 


(E\0) 


1 


2 


3 


(E\0) 


-2 





4 


(G 4 |r) 


t-ri \r\\ 


U 


(G 4 |U) 


1 
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n 

z 


1 
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(G 4 |UJ, (G 4 |UJ 
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o 
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(C 2 |r) 


(tCa|0) 
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(£|0) 


1 





1 


(E\0) 


-2 
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(i|T) 


(B|0) 


6 


(E\Q) 


1 





3 


(E\0) 


2 





-6 


(iC||T) 


(C2|0) 


2 


(£|0) 


1 


4 


3 


(E\0) 


-2 
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(iG 4 |0) 


(C 4 |r) 


2 


(G 4 2 |0) 
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(G 4 2 |0),(G 4 2 |0) 











(tCa|0) 


(C 2 |r) 





(£|0) 
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2 


1 


(E\0) 


-2 









Only classes that have non-trivial characters are shown. 
Ai ® Ai = Tf © rjj e and A 6 ® A 6 = rf e rf 2 e 2r^ 5 e r^ 5 without time-reversal symmetry 



TABLE X. Characters of the product representations at A = 2ko. 
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Xa 18 ' 


\l Ai®Ai 
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Q is the identity element. A6 ® A6 = Ai © A'i © As without time-reversal symmetry 



TABLE XL Characters of the product representations at E = fco(l, 1,0). 
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N(QC) 
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Ai ® Au = Si © E 4 and A 6 ® A 6 t = 2Ei © 2E 2 © 2E 3 © 2E 4 without time-reversal symmetry 



TABLE XII. Characters of the product representations at L. 
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Q is the identity element. Only classes that have non-trivial characters are shown. 
L t ® L~l = F~t © © r+ 2 © 2r^ 5 © 2rJ 5 without time-reversal symmetry 



TABLE XIII. Characters of the product representations at X. 
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Four other non-trivial classes that have exactly the opposite characters are not shown, 
ii ® Lit = Xi ® Xs and ® L^ t = 2Ai © 2A 2 © 2A 3 © 2A 4 without time-reversal symmetry 



